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ABSTRACT 

We have found relations between the temporal and spectral properties of 
radiation Comptonized in an extended atmosphere associated with compact 
accreting sources. We demonstrate that the fluctuation power spectrum density 
(PSD) imposes constraints on the atmosphere scale and profile. Furthermore, 
we indicate that the slope and low frequency break of the PSD are related to 
the Thomson depth r of the atmosphere and the radius of its physical size 
respectively. Since the energy spectrum of the escaping radiation depends also 
on t (and the electron temperature kT e ), the relation between spectral and 
temporal properties follows. This relation allows for the first time an estimate 
of the accreting matter Thomson depth r independent of arguments involving 
Comptonization. We present figures for the light curves and PSD of different 
energy bands, the photon energy spectra and the phase lags as functions of the 
variability frequency. The temporal properties of the high (soft) and low (hard) 
state of black hole sources are discussed in this context. 

Subject headings: accretion — black hole physics — radiation mechanisms: 
Compton and inverse Compton — stars: neutron — X-rays: QPO 
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1. Introduction 

Much of the information about the conditions in the vicinity of compact accreting 
objects, the region where most of the accretion kinetic energy is themalized and converted 
into radiation, has so far come from observations and modeling of their energy spectra. 
Specifically, the observations of power-law high energy (E > 2 — 50 keV) spectra has 
established the process of unsaturated Comptonization of soft photons by hot electrons 
as the main process by which these spectra are produced (see e.g. Sunyaev & Titarchuk 
1980). However, the spectral information alone does not suffice to determine the physical 
parameters of the system, since it can only determine the combination of the Thomson 
depth, To, and the electron temperature, T e , which is relevant to the Comptonization 
of soft photons by the hot electrons, namely the Comptonization parameter y defined 
as y ~ (AkT e /m e c 2 )TQ, for r 3> 1. In this respect, one should bear in mind that for 
quasi-spherical, free-fall accretion, the above parameters, i.e. r , y, are independent of the 
mass of the accreting object, provided that the accretion proceeds at the same fraction of 
its Eddington value. Thus, in order to completely determine the physical parameters of 
these systems, additional, independent information is needed to set the scale of the radius 
and the density of the emitting region, information which can only come from the dynamics 
of accretion. These arguments have motivated the study of time variability of accreting 
sources and in particular of galactic black hole candidates; it was thought that the absence 
of a solid boundary in these bright sources would allow one to probe the dynamics of 
accretion onto the black hole. 

However, the study of time variability of black hole candidate sources, in particular 
the archetypal galactic source Cyg X-l, has yielded a number of results which were rather 
unexpected when viewed within the standard framework of the viscous accretion disk 
model for the dynamics of matter accretion and emission of radiation from these sources: 
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(a) The fluctuation power spectral densities (hereafter PSD) are generally power laws of 
indices s ~ 1 — 1.5 in the variability frequency uj, i.e., |F(u;)| 2 oc uj~ s for uj > uj c and 
|F(u;)| 2 oc uj° for uj < uj c (see e.g. van der Klis 1995). The values of the slope, s, of the PSD 
are thus significantly flatter than these expected from exponential shots associated with the 
turbulent dynamics of the accreting gas near its last stable orbit (s > 2). Furthermore, 
the turnover frequency, uj c , is generally several orders of magnitude smaller than those 
associated with the dynamics responsible for the emission of high energy radiation. For 
Cyg X-l, the archetypal galactic black hole candidate, s ~ 1 (flicker noise), while uj c ~ 0.05 
Hz, a far cry from the kHz frequencies expected on the basis of dynamical considerations. 

(b) In addition to these broad band spectral features, the PSD of black hole candidates 
also exhibit QPO's. The observation of these features, thought to result from the "beating" 
between the Keplerian frequency of an accretion disk and the rotation of an underlying 
neutron star, poses a problem for systems thought to contain black holes, (c) The spectral - 
temporal studies of Miyamoto et al. (1988) of the time lags between hard and soft photons 
as a function of the variability frequency uj, has indicated that these lags, at least in Cyg 
X-l, decrease with increasing variability frequency uj, thus ruling out Comptonization by 
an electron cloud of uniform density and temperature. 

Motivated by the discrepancy of the above systematics of time variability in black hole 
candidates from those expected on the basis of the "standard" dynamical models of viscous 
accretion disks, we present in this study an alternative model for the time dependent 
spectral formation of Comptonized radiation which can reproduce the basic observational 
features described above. In §2 we describe our model and demonstrate the relation between 
the PSD shapes and the slopes of the photon energy spectra; we also emphasize the nature 
of the white noise spectra below uj c as well as the physics associated with the low value of 
this turnover frequency. In §3 we discuss the issue of the very low-frequency noise and we 
briefly review the relevance of the present model to RXTE observations of Cyg X-l in its 
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low and high states. Finally, we discuss and summarize our results in §4. 

2. The Extended Atmosphere Model 

Our study uses as background the model of Chakrabarti & Titarchuk (1995; 
hereafter CT95), i.e. quasi-spherical accretion with a centrifugally supported shock, which 
thermalizes the accretion kinetic energy and gives rise to the hot electrons responsible 
for the Comptonization of soft photons. However, we postulate, in addition, that the 
quasi-spherically accreting component behaves like a hot "atmosphere" of constant electron 
temperature T e ~ 50 keV and, more importantly, with a density profile n(r) oc 1/r in 
radius, extending to r = r c ~ 10 4 R s , where R s = 2GM/c 2 is Schwarzschild radius. This 
atmosphere is thought to be the result of preheating of the accretion flow by the radiation 
produced at the shock, as discussed by Zeldovich & Shakura (1969) and Chang & Ostriker 
(1985 and references therein) (see additional discussion on the issue of this specific density 
profile in §4). 

We assume, for simplicity, that the geometry of this configuration is spherical and that 
there is a source of soft photons within the spherical shock boundary of radius r s h. The 
electron density is considered to be constant, n+, inside this boundary and n + r sh /4r in the 
atmosphere outside. The discontinuity of a factor 4 in density across the boundary r sh is to 
account for the density jump across the shock. The physical size of this cloud is determined 
by the total optical depth r . We have calculated the response of this configuration, i.e. 
the region interior to the shock plus the extended atmosphere, to an impulsive input of soft 
photons within the radius of the shock. The calculations were carried out by a modified 
Monte Carlo code based on the method described in Hua & Titarchuk (1995). Figures la 
and lb show the resulting light curves at different energy bands for clouds with r = 1, 2 and 
3, n + = 1.6 x 10 17 cm -3 , r s h = t x 10~ 4 light seconds and outer radius of the atmosphere 
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r c = 5 x 10 3 r sh ~ r 0.5 light seconds. It is apparent that their shapes have the form of 
power laws over the time range 1CT 3 s to ~ 1 s, followed by an exponential cutoff at times 
of order a few seconds. For comparison, the exponential light curves from a cloud with 
uniform density n = 2 x 10 14 cm -3 , radius r = 1.5 x 10 10 cm (~ 0.5 light seconds) the same 
electron temperature T e and r = 2 are also shown (dotted curves). 

The power law form of the light curves, /(£), is the result of photons scattering in the 
extended atmosphere whose optical path has a logarithmic radial dependence 

rr 

r= a T n e dr = r sh [l + 0.25 • ln(r/r sh )] for r > r sh (1) 
Jo 

where the electron density n e = n + for r < r sh and n + r sh /4r for r > r sh ; ar and r sh are the 
Thomson cross section and optical thickness of the shock respectively. In the case that the 
optical depth of the extended atmosphere is a few, most photons escape either unscattered 
or after one scattering (we exclude in these considerations the photon spread in time over 
the short time scales ~ r^/c associated with the uniform sphere at the center of radius 
r s h). The light curves of figure 1 present the distribution in the delay times of the escaping 
photons with respect to the light crossing time, for photons which have scattered once. The 
light curve is related to the distribution of photons over the time of the first scattering t\. 
The probability of scattering in an interval between t\ and ti + dt\ is P{ti)dti = /(r)dr 
where /(r) is the probability density to scatter between depths r and r + dr. In our case 
this probability density is constant and equal to l/r , i.e. inversely proportional to the 
total depth of the atmosphere. Therefore, 

. , 1 dr 1 dr dr 1 . . 

P(ti) = --rr- = --TTT = -<™ e (r)c (2) 
r dti r dr dti t 

For a density distribution which is a power law, i.e. n e (r) oc r~ a (we herein consider 
the case a = 1) P(ti) oc One can establish, by looking at the geometry of photon 
trajectories in a spherical geometry of radius R, that for a photon which has scattered at 
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radius r x the average distance before escape is given by 



= o r i 



l Ti? l /i? 2 




indicating that for small values of T\ the average delay time td of photon escape after one 
scattering relative to the light crossing time of the entire atmosphere of size R is 



The linearity between the scattering time t\ and the delay time td allows us to substitute 
td in equation (2) above, leading to a light curve which is a power law in t (or better in td) 
with index equal to that of the density profile index a. 

The Monte Carlo calculations shown in figure 1 support the arguments presented 
above. The slope of the the light curve for the uniform sphere of small scattering depth is 
close to zero (dotted curve) and that of the atmosphere considered here is close to -1; we 
have also used a density distribution appropriate to free-fall i.e. n(r) oc r~ 3 / 2 with results 
in accordance with the above arguments. It is well known (Sunyaev & Titarchuk 1985) that 
for photons escaping from a source of finite size, the distribution of number of scatterings is 
an exponential i.e. f(u) oc exp(— f3u), where u is the scattering number and the parameter 
/3 is inversely proportional to the mean number of scatterings, u, which is greater than 1 
for clouds with total optical depth r >l. Thus, as expected on the basis of this argument, 
the light curve has an exponential turnover at delays longer than 1/(5 times the scattering 
time at the largest decade of radius. These arguments therefore present the possibility that 
measuring these light curves, e.g. by computing the resulting PSDs, can provide a tool for 
uncovering the density profile of the atmosphere . 

As the optical depth of the atmosphere increases, multiple scattering (diffusion) effects 
in the extended atmosphere become more important and the above arguments, which rely 
on the single scattering approximation, have to be reformulated. The power law indices of 




(3) 
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the light curves evolve to smaller values; in fact, in the uniform sphere case and for very 
large optical depths, as well known, the light curves become rising in time (see e.g. Hua & 
Titarchuk 1996; hereafter HT96). The flattening of the light curves is indicated clearly in 
our our Monte Carlo calculations presented in figures la and lb. This change in the light 
curve slopes reflects directly on the resulting PSD spectra, shown in figures lc and Id. It 
can be seen there that the PSD curves are also power laws in Fourier frequency u, since the 
Fourier transformation of a power-law, t~ a , is also a power-law in frequency uj a ~ l . Thus the 
flattening of the light curves leads to steepening of the PSD spectra, which however, as a 
rule, are flatter than those corresponding to the PSD of a uniform source (oc u~ 2 ; dotted 
curves), which correspond to exponential-shot light curves. 

Assuming that the temperature of the extended "atmosphere" is constant, a change in 
its total optical depth r , which results in a change of the light curve index and a steeper 
PSD, will manifest itself, in addition, in the energy spectrum of the escaping radiation: The 
increase in optical depth will lead to a harder energy spectrum for the emerging photons. 
Therefore, our considerations, supplemented by the Monte Carlo calculations, support 
to the conclusion that, under the assumption of the existence of an "atmosphere" with 
the density profile described above, there should exist well defined correlations between 
the slopes of the escaping photon spectra and those of the PSD. We believe that RXTE is 
uniquely suited to search and document the presence of such correlations. 

In Figure 2 we plot the energy spectrum (solid curve) resulting from the extended 
atmosphere with parameters kT e = 50 keV and r = 3. It is seen that it is almost identical 
to that from a uniform plasma cloud with the same electron temperature T e but a different 
optical thickness r = 2 (dashed curve). For comparison, we also plot the energy spectrum 
from a uniform cloud with the same temperature and r = 3 (dotted curve). It is seen 
that this spectrum is harder than the other two. This is because for the same total optical 
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thickness (r = 3), photons in the uniform cloud find it harder to escape than those in 
the 1/r atmosphere. This example shows clearly that spectroscopic analysis alone can not 
provide complete information about the source structure and to get a complete picture of 
the spectral formation it is necessary to analyze in conjunction the temporal as well as the 
spectral data. 

The time information is of course of vital importance for setting the physical scale of the 
system. Since it is believed that Comptonization is the main mechanism for the production 
of X-rays, measurements of the time lags between two different energies in the X-ray band 
should give us a direct estimate of the density of the region at which Comptonization takes 
place. Based on simple arguments concerning the dynamics of accretion onto compact 
objects, the associated time lags are expected be of the order of msec for galactic objects. It 
is worth noting that, for the uniform density hot clouds considered almost exlusively todate 
in the literature, the time lags between soft and hard photons are constant and independent 
of the Fourier frequency uj. Extensive caculations of this type (i.e. of scattering in clouds 
of uniform density) are given in HT96. For the non-uniform density profile of the extended 
atmosphere discussed herein, the shape of the phase lag as a function of the variability 
frequency uj also changes. In Figure 3, we plot the phase lags of photons in the 10 - 20 
keV band relative to those in 1 - 10 keV band, for the four cases depicted in Figure 1, as 
a function of a;. As in Figure 1, the dotted curve corresponds to the configuration of the 
uniform cloud described earlier and its shape is similar to those displayed in Figures 14 - 17 
of HT96, except that here the curve extends to higher frequencies with a roughly constant 
slope instead of a sharp drop. This is because in present analysis we use much finer time 
bins than before. However, the maximum at frequency ~ n e a T c/T (HT96, Eq. 11) is 
clear and sharp. On the other hand, the phase lags corresponding to the 1/r atmosphere 
configuration (solid curves) are much flatter and there is no maximum present across the 
entire range of frequencies. It is interesting to point out that this type of phase lag is hinted 
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in the GINGA observations of Cyg X-l (Miyamoto et al. 1988) and also seen in the recent 
observations of the high state of Cyg X-l (Swank et al. 1996). 

As pointed out in HT96, the PSD associated Comptonized photons of a given energy, 
depends strongly on the energy of the source, soft photons. In Figure 1, the source photons 
have a blackbody temperature 2 keV. In Figure 4, we present the PSD curve (solid curve) 
in the energy range 10 - 20 keV resulting from the same configuration as that of r = 2 in 
Figure 1, but with the source photons at blackbody temperature kT = 0.5 keV. Compared 
to the corresponding PSD curve in Figure 1, it is seen that PSD with lower source photon 
energy is steeper, in agreement with the analysis of HT96. 

In Figure 4, we also present one more PSD (dotted) curve, corresponding to the same 
light curve as the solid one but for different time ranges and time bins. For the solid curve 
the light curve is calculated over a range of 4 seconds in 4096 bins. For the dotted curve, 
the light curve is calculated over 32 seconds in the same number of bins. It is seen that the 
PSD turns flat for frequencies below u c ~ 0.25 Hz while keeps parallel to the solid one above 
uj c except at the highest frequency end, which is obviously due to aliasing (see e.g. Press et 
al. 1992). It is found from examining the light curve that the time scale corresponding to 
l/u c = 4 seconds is actually the time scale of the light curve beyond which the latter drops 
virtually to zero. Thus we have found a possible physical meaning for the shoulder (break) 
frequency u c in the PSD curves, which is common in the PSD of many sources, namely, 
it indicates the time scale of the light curve, or alternatively the size of the extended 1/r 
atmosphere. The "white noise" below u c reflects the average frequency of the shots while 
the power-law above ui c reflects the time structure within one single shot. 
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3. The very low frequency noise 



In addition to the fluctuation PSD described above, galactic X-ray sources exhibit 
power-law type PSD at frequencies much lower than u c , so that the white noise component 
at uj < uj c appears only as a shoulder in the much broader PSD spectrum. We have thus 
reasons to believe that this very low frequency component in the PSD is separate from 
the one described above. One should bear in mind that at radii r ^> r c , related to the 
dynamics associated with this component in the PSD, the effects of preheating are thought 
to be small and accretion proceeds in the form of the standard viscous a-disk (Shakura 
& Sunyaev 1973 hereafter SS73). As a result, one can consider the fluctuations in the 
accretion rate M = M + 5m to be resulting from small random variations in the viscosity 
parameter a, i.e. a = a + 5. This model was considered by Lyubarskij (1995), who, 
by solving the angular momentum diffusion equation, showed that the resulting accretion 
rate and hence the luminosity, has a flicker-type fluctuation spectrum, provided that the 
characteristic time of the fluctuations in a are of the order of the viscous time scales, and 
the fluctuations in a at different radii are uncorrected. 

We provide here a heuristic derivation of this fact for the sake of completeness of the 
discussion. The accretion time r a from a radius r in an a-disk is 



The variation 5 in the viscosity parameter a in an annulus leads to the mass accretion 
rate variations 5rhi oc 5(r, t). In fact the contribution in the mass variation caused by 
the viscosity variations is proportional to the annulus area rdr, the surface density u and 
inversely proportional to the accretion time r a i.e.. 




2 



-i -i 



(4) 



Sihi oc 



5 ■ urdr 



(5) 



a 
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In different accretion disk regions (SS73) (region a: p r ^> p g and a T ^> a - //, region b: 
Pr <S p g and o"t 3> cr//, region c: p r <C p 9 and er T <C o - //) there are different dependences of 
the surface density u and thickness of disk H on radius. For region a, u oc r L5 , i7 is almost 
constant and hence r a oc r 3 - 5 (see Eq. 4) ; for region b, u oc r~ a6 , iJ oc r 21 / 20 and r a oc r L4 ; 
for region c, u oc r 75 , oc r 9 / 8 and r a oc r 3 5 . Despite the differences in the behavior of the 
main parameters over the disk one can check by substituting the above expressions of u, H 
and r a in Eqs. 4 and 5 that the mass (or luminosity) variations are the same for all zones 
a-c, namely 

at a- dr dTa duJ (R\ 

dL oc ami oc — oc — = — (o) 

r r a u 

and thus the resulting power spectrum p{u) is 

oc u~ l . (7) 

In other words if the amplitude of the variations in a are the the same at different radii 
then the amplitudes of accretion rate (luminosities) variations are the same at different 
time scales (Lyubarskij 1995). 



4. Discussion and Conclusion 

We have presented above a model which implies corelations between the spectral and 
temporal properties of accreting compact objects and have pointed out that this relation, 
along with independent measurements of the electron temperature T e , can allow a complete 
specification of the physical parameters of these systems, including the density profile of 
the extended atmosphere. 

Much of the present discussion relies on the specific form of the density n(r) of the 
accreting matter as a function of radius r. Clearly, this is not the free-falling solution 
customarily used in association with spherically symmetric accretion, so a few comments 
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are in order. If the presence of the extended atmosphere is due to the effects of preheating, 
as suggested in §1, then at the edge of this atmosphere, at r ~ r c , one would expect the 
random and rotational velocities of matter to be comparable and that the subsequent 
evolution of the accreting fluid in radius to be predominantly governed by the removal 
of its angular momentum. The agent responsible for this process is considered herein to 
be the interaction of the fluid with the photons produced at the base of the extended 
atmosphere, i.e. near the Schwarzschild radius R s . In order to estimate the effectiveness of 
this process we compare the free-fall time scale, tff, with that of the viscous time scale, t V i SC , 
assuming that the density has the required profile. Assuming that the density profile of the 
atmosphere has the form n(r) = n (r /r) between radii r and r c , with r c ~ 10 3 — 10 4 r , 
mass conservation dictates that the the infall velocity, v(r), will also have a similar scaling 
between tq and r c , i.e. v(r) = vo(r /r), yielding for the free-fall time tff ~ r/v{r) = r 2 /vqVq. 
The photons can achieve the required angular momentum removal if tff ~ t V i SC . 

Under the assumed density and velocity scalings, the angular momentum of the fluid 
at a radius r would be C ~ m p n(r)v(r)r ■ r 3 . This fluid interacts and transfers its angular 
momentum to photons diffusing from the interior regions of the atmosphere which carry 
substantially smaller angular momentum. The specific density profile we have assumed, 
which implies an equal optical depth per decade of radius, guarantees that: (a) The angular 
momentum of the photons, as they traverse each decade in radius, is much smaller than 
the local angular momentum of the fluid, (b) The photons interact with the fluid at every 
decade in radius with probability ~ 1 before their escape, carrying away angular momentum 
at a rate C ~ (L/c 2 )v(r)r, where L is the total luminosity of accretion. The viscous time 
scale for the fluid at radius r is then the time required for the fluid at radius r to get rid of 
its angular momentum, i.e., 
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C n(r)c 2 r 3 m p . . 

tvisc — — \PJ 

The luminosity L is given simply by the total accretion rate M and the innermost 
radius of the accretion disk r\, i.e. L ~ (3Mc 2 {R s /rj) ~ j3m p noVor 2 c 2 (R s /ri) with the 
parameter f3 allowing for accretion onto the compact object other than that associated with 
the spherically accreting component (most likely associated with a viscous geometrically 
thin disk). Substituting the expression for L onto equation (8) above yields for the viscous 
time scale 

r 2 1 1 

The above expression indicates that the removal of the angular momentum of the accreting 
matter by the photons produced near r ~ fi s < r c can in fact proceed on dynamical time 
scales and thus preserve the assumed density profile provided that {R s /ri)(5 ~ 1, i.e. that 
the total accretion rate (including an additional component, most likely from an accretion 
disk) is ~ (ri/R s ) times that associated with the spherically symmetric component. We 
reiterate that, most likely, this process is possible only for the density profile prescribed 
above (n(r) oc 1/r), since this is the only profile which allows for significant photon 
scattering, and hence removal of angular momentum, from a large range of radii. 

The presence of the extended atmosphere discussed above, due presumably to upstream 
heating of the accreting matter, may bear relevance not only to accreting black holes as 
discussed herein, but also to accreting neutron stars, with the overall similarity of their PSD 
spectra (Van der Klis 1995) not being simply coincidental. It is interesting to speculate that 
coherent oscillations of this configuration (the extented atmosphere) may in fact be related 
to the QPO phenomenon, a phenomenon prevalent among all members of the more general 
class of accreting compact sources, whether neutron stars or black holes. The observed 
QPO frequencies are generally smaller than those associated with the dynamical motions 
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near a black hole or a neutron star, and similar to those associated with the outer radius 
r c of our model, implying the possible relevance of this characteristic time scale in these 
systems too. This model may also bear relevance to the recently discovered kHz QPO's in 
four LMXBs (Strohmayer et al. 1996, Zhang et al. 1996), as this is the time scale expected 
for variability in the centrifugally supported shock presumably present at the base of the 
extended atmosphere, especially in objects containing neutron stars rather than black holes 
(Titarchuk k Lapidus 1996). 

The presence of a reflection-by-cold-matter component and fluorescent K a line emission 
in the observed spectra of black hole candidates and AGNs (e.g. Ebisawa et al 1996) could 
also be related with an extended atmosphere of the type discussed above: A significant 
fraction of hard photons can in the present case random walk over the entire extended 
atmosphere and the possibility arises that it can be reflected by far-away, relatively cold 
parts of the visous accretion disk to produce the required components. 

In addition to the above and in accordance with the model of CT95, one would expect 
a variation in the temporal properties of accreting compact sources while at their different 
spectral states. As discussed by CT95, the changes in the spectral states derive from the 
presence of a sufficiently large number of soft photons from the viscous accretion disk to 
cool the electrons to temperatures at which Comptonization is ineffective. This manifests 
as a prominent thermal-like peak at energies ~ 1 keV. Under the same conditions, one 
could expect the extended hot atmosphere to be absent, a fact which should manifest itself 
in the associated PSD by the absence of power at the lowest frequencies. 

We have presented above a model which implies the existence of a correlation between 
the temporal and spectral properties of accreting compact sources. Our model is consistent 
with certain general characteristics associated with the PSD of accreting compact sources 
observed to date, in particular with their spectral forms and the breaks associated with the 
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low frequency turn-overs. It is also consistent with the observed frequency dependent phase 
lags between photons of different energies and can also account for their photon spectra. 
This model is specific enough that we believe can it be meaningfully tested by combined 
spectal - temporal measurements of these sources. 
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Fig. 1. — The light curves and PSD for the extended atmosphere of temperature kT e = 50 
keV described in the text. Three cases of total Thomson optical depths r = 1,2 and 3 are 
shown. The electron density has the form 0.25n + r s h/r for radius r > r s h and n + for r < r s h, 
where n + = 1.6 x 10 17 cm' 3 and r sh = t x 10~ 4 light seconds. For comparison, the cloud 
with uniform density n = 2 x 10 14 cm -3 with the same temperature and r = 2 is also shown 
(dotted curves). The left two panels show the light curves and PSD for emissions in the 
energy range 1 10 keV, while the right two panels for 10 20 keV. 

Fig. 2. — Comparison of the energy spectra of emissions from uniform clouds and extended 
atmospheres with the same temperature and optical depth r = 3. The solid curves 
indicate the energy spectrum resulted from same extended atmospheres as in Figure 1. The 
dotted curves indicate those from uniform clouds. The spectrum resulting from an extended 
atmosphere of r = 3 (solid curve) is almost indistinguishable from that of a uniform cloud 
with r = 2 (dashed curve). Both spectra have the same parameter (3 = 0.334 although they 
have different r . 

Fig. 3. — The hard X-ray phase lags resulting from the same atmospheres as in Figure 1. 
The dotted curve indicates the phase lag resulted from a cloud with tq = 2 and uniform 
density n = 2 x 10 14 cm -3 . 

Fig. 4. — The power spectra resulting from the same shot light curve but different time 
ranges. The emissions are from an extended atmosphere with r = 2 and kT e = 50 keV. 
The density profile is similar to the one with r = 2 in Figure 1, but with source photons at 
blackbody temperature 0.5 keV. The solid curve represents a time range of 4 seconds while 
the dotted curve for 32 seconds. 
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